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I. INTRODUCTION 

At low temperature, the magnetic response of small rings has a component which changes 
periodically with the applied magnetic flux |l|-^. Although this phenomenon, i.e. the ex- 
istence of a persistent current, was predicted many years ago 0-0], the magnitude of the 
effect is not well understood: The experimentally determined persistent currents in metallic 
rings [|I]-0] are much larger than theory P|-PU|. predicts. For systems where the electron 
motion is diffusive, i.e. the mean free path I is much smaller than the circumference of the 
ring L, important theoretical results for the average persistent current are the following: For 
non- interacting electrons, the current is of the order / ~ A/0O) where A is the mean level 
spacing at the Fermi energy, A = (2A/oV)^^, with A/q the density of states per spin and V 
the volume, and (po = h/e is the flux quantum |TT|-[T3|. This result applies for temperatures 
below the Thouless energy, which is given by Ec = hD/L'^, D = vpl/3, but large compared 
to the mean level spacing. For interacting electrons, the Coulomb interaction contributes to 
the persistent current [|^,0, namely / ~ ^*Ec/(f)o, where fi* characterizes the strength of 



the interaction. For the metallic rings in the experiments the inequalities A -C kT, Ec were 
fulfilled. 

The experimental results can be fitted almost perfectly with the 'interacting' theory of 
| I^ if we choose /i* = ±0.3. However the precise value of the interaction constant is not 



known from theory. To first order in the screened Coulomb interaction, the interaction 
constant is of the correct order compared to experiment, however, a repulsive interaction 
scales downwards when higher order corrections are taken into account p7| , p8i [T5|. For 



example, for copper, the material of which the average persistent current has been measured. 



/i* has been estimated |]T5[ to be smaller than the bare value of 0.3 by a factor of about five 
to ten. 

Many recent studies of persistent currents were devoted to the understanding of the 



interplay of disorder and electron-electron interaction [p!^ -p5 



In the present article we study a new mechanism which may induce a persistent current. 



namely magnetic impurities. This is a novel phenomenon: While the sensitivity of the 



persistent current to magnetic scattering has been predicted [p3| , |29| -[3T[] in previous studies, 
we point out that these investigations started from the theory of persistent currents without 
magnetic impurities, and introduced magnetic defects as an additional perturbation. In such 
an approach, an enhancement of the persistent current is never to be found. 

In contrast, in our recent article where we concentrated on persistent current fluctuations, 
we already concluded that there should exists a new contribution to the average persistent 
current induced by magnetic impurities. Let us recall the argument |3l| : In a diffusive ring, 
the stochastic fluctuations of the persistent current are of the order ((5/)^) ~ {Ed (poY- 
Changes in the impurity potential lead to a variation of the persistent current with this 
order of magnitude. Considering magnetic impurities as a spin-dependent scattering po- 
tential, which changes slowly due to spin relaxation processes, one finds temporal cur- 
rent fluctuations of the order (((5J)^) ~ {Ec/ipoY o^ ^^e time scale of the spin relaxation 
time tk- Using the fluctuation-dissipation theorem we relate these current fluctuations to 
the dynamic linear response x{'^) to a time-dependent magnetic flux 0(u;) with the result 
Imx(co' ~ I/tr-) ~ {Ec/ (poY {h/ Tk) / kT . This suggests that at low frequency, there is also a 
contribution to the real part of x of the order x(u; — *> 0) ~ {Ec/(po)^/kT. In the static limit 
X is the derivative of the persistent current with respect to the flux: x{^ = 0) = d^I{(f)). In 
the present paper, we conflrm this earlier suggestion ||3l[] explicitly, namely / ~ El/{kT(f)Q). 

In section II, we start with a few simple considerations, which demonstrate that a flux- 
dependence of the free energy induced by magnetic impurities may indeed exist. Then we 
extend the calculations and discuss (1) the role of spin-orbit scattering, (2) strong spin- 
flip scattering due to high concentrations of magnetic impurities, and (3) strong spin-flip 
scattering due the Kondo effect. 



II. WEAK SPIN SCATTERING 

We start from the free energy of a spin free spin [S = 1/2), given by 

n = -kT\n (2 co^h^^\ (1) 

The spin is coupled to the conduction electrons; we assume the usual local exchange Hamil- 
tonian H = Js(R) ■ S, where s(R) is the spin density of the conduction electrons at the 
impurity site. To first order in J, this coupling is equivalent to an additional magnetic field 
for the impurity spin, and thus leads in Eq. (|l|) to the replacement 2jibH — > 2jibH + Js^{^). 
For electrons on a ring, the local spin density depends on the magnetic flux cf) penetrating the 
ring. As a consequence the free energy is flux-dependent, and a persistent current / = —d^Vt 
exists. 

For the explicit calculation, we decompose the electron spin density in its fiux- 
independent mean value plus the flux- dependent fluctuations, s^(R) = (s^(R)) +(5s^(R, 0). 
Then we expand VL with respect to 5s^ and average over disorder, with the result 
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(fi) = fi((s^)) - ^x'^.f{5s\K)5s\K)) + ■ ■ ■ , (2) 

where x^^ = —djjQ is the (longitudinal) susceptibility of the impurity. These rather simple 
considerations already lead to a novel contribution to the persistent current. 

Equation (H) is valid for an arbitrary geometry. We consider a ring of circumference L 
and transverse dimension L± <^ L. We assume diffusive motion of the electrons, i.e. the 
(elastic) mean free path is much smaller than the circumference of the ring, / <^ L. 

We evaluate the fluctuations of the local spin-density using the standard Green's func- 
tions technique. We average over impurity configurations, keeping only the diagrams with 
one particle-particle ladder, i.e. one cooperon, as shown in Fig. |T[ The flux sensitivity is 
determined by the long-wavelength, low-frequency contribution of the cooperon. The result 
is{h=k = l) 

^"^^ e,tJ k,k',q s,s' 

xG,(ie, k)GAK(^ -uj),-k + q)a(ie, k')GAK^ - ^), -k' + q), (3) 



where V is the volume and Gg the averaged Green's function for electrons with spin s, and e 
and uj are odd and even Matsubara frequencies. The relevant components of the cooperon, 
including the Zeeman effect, are given by 

C++++ = C____ = ^-^^^^—^ (4) 

provided e(e — cu) < 0; otherwise no singular contribution arises. As usual, Us = ^fisH 
denotes the Zeeman splitting; see also App. A, where we summarize the results for the 
cooperon including spin-flip and spin-orbit scattering. For simplicity, we consider the limit 
1/ts = 1/tso = for the moment. The diffusion constant, D, is given by Z^ = Vpr/S, where 
T denotes the elastic scattering time. 

Equation (j^) is similar in form to the expression for the fluctuations of the local electron 
density, which have been evaluated in |T^ in order to determine the Coulomb interaction 



contribution to the persistent current. In ||T^ both Hartree and exchange contribution to 



the persistent current are calculated. The results presented here and below correspond to 
the Hartree terms; we neglect the exchange contributions, which are insignificant in most 



cases. In order to evaluate Eq. @ we rely on the results given in [|l4|. 

For the ring geometry, the transverse dimensions of q are frozen out, such that it can be 
considered one-dimensional, assuming the values g = 27r(n + 2(f)/(j)o)/L, where 0o = h/e is 
the flux quantum and n is an integer number. The summations over k and k' are converted 
into integrals, and from each integration we find a factor (27rA/or). The e summation is then 
feasible since the e dependence is only due to the condition that e(e — uj) < 0. Finally we 
arrive at 

((5s^(R)(5s^(R)) = ^T V ^ V ( ^-— - Re ^\ . ] (6) 

The spin fluctuations are periodic in with periodicity h/2e, and can thus be represented 
as a Fourier series, 

(5s^(R)5s^(R)) = Y. ^mexp(27rim0/0o). (7) 

m=—oo 



Note that only the even components are non-zero, with the result 



V ^0 7-00 \uj + 4ti^Dx^/L^ uj + ATiWxyL^ + iuJsJ 

(8) 

All even coefficients are related to the second coefficient at rescaled temperature and Zeeman 
energy, A2m(T,Us) = A2{m'^T,m'^u!s)/\m\^, so it is sufficient to calculate A2. Since Am = 
A-rn only terms proportional to cos(27rm0/0o) survive in Eq. (J^). The Fourier series of the 
grand potential is 

00 
(Q) = {Qo) + E (^-) cos(27rm0/0o), (9) 

m=l 

with (Qm) = ~X^^J^Am/4:fi%. Finally, the Fourier components of the persistent current are 
given by (/„) = em{Qrn)- 

Numerical results for the ffist non- vanishing Fourier component of the grand potential, 
(^^2), as a function of temperature and for different values of the Zeeman energy are shown 
in Fig. 1^. For a finite but small concentration of magnetic impurities, we have to multiply 
the result with the number of impurities A^^- In Fig- § we choose the concentration such that 
the spin-flip scattering rate equals the Thouless energy: l/r^ = 2-nnsN'Q.PS{S + 1) = E^. 
We find that (^2) goes to zero both in the limit of a weak and a strong magnetic field: For a 
weak magnetic field, there is no electron spin polarization and consequently no fluctuations 
in the spin polarization, while for a strong magnetic field {ug ^ T), (0.2) goes to zero since 
the susceptibility goes to zero. The maximum value is found for intermediate values of the 
magnetic field, where there are fluctuations in the electron spin density but the impurity 
spins are not yet fully polarized, such that x"^ ¥" 0- 

III. STRONG SPIN SCATTERING 

In the case of strong spin-orbit scattering or strong spin-flip scattering, there are fluctua- 
tions in the local electron spin density even without Zeeman effects and we find a significantly 
enhanced persistent current. 
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A. Strong spin-orbit scattering 

Let us consider spin-orbit scattering, assuming that spin-flip scattering is weak. In this 
case it is straightforward to generalize the calculations described above. We take into account 
that there exist fluctuations in the electron spin density not only in z-direction, but also in 
X- and y-direction. Equation (^ generalizes to 

(fi) = n{{s^)) - ^.f ir^{5s^5s^) + Ax'-iSs-Ss^)) , (10) 

with the susceptibilities given by 

' -X- = 1^ f 1 - tanh^ ^) (11) 



4/^1'" AT V 2T. 

-^^X = tanh -^. (12) 

The fluctuations in the local spin density are 

(^Ss^Ss^) = ^ E ^(2^-^or)^ E <^:'s'Css'ss' (13) 

w>0 g,s,s' 

i^^^S^') = V ^ f (27rAror)2E^-++^- (14) 

with the expressions for the cooperon given in App. A. 

In the limit of zero spin-orbit scattering, C_++_ is zero and thus the transverse spin 
fluctuations are zero. Then the longitudinal fluctuations are identical to the expression 
given in Eq. (P). For very strong spin-orbit scattering, on the other hand, (I/tsq ^ E^ and 
cuf <^ Ec/tso), the cooperon simplifies considerably: C++++ ~ 0, 27rA/or^C+_+_ ~ l/2A^o, 

and 27rA/or^C-| ^ ~ — l/2A^o, where A^o = |^| + Dq^. The flux- dependent part of the 

thermodynamic potential is then given by 
AT, 



8/ii 



J2 (^x""{Ss'6s') + 4x+"(5s-(5s+)) 



2 
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^J.u + Dq- ^''^ 



Fig. ^ shows the amplitude of the second harmonic, (^2), for the same values of ujs, T, 
and 1/ts as in Fig. |^. Spin-orbit scattering changes the sign of {^2), it does not go to zero 



for cus — ^ 0, and especially the amplitude of (^^2) is much larger than without spin-orbit 
scattering. For low temperature, the fiux-dependence due to magnetic impurities dominates 
the fiux-dependence due to the Coulomb interaction. The latter for low temperature is given 
by |T^ (^22) = (4/7r)yU*i?cexp(— T/3i?c)- In the figure we use the estimated value for the 
effective interaction, /i* ^ 0.3/5 = 0.06. 

We note that for the temperatures in the figure, the fiux-dependence of the grand po- 
tential due to magnetic impurities is well approximated by 

(^2) = -^ fl -tanh^^ + — tanh^) 3_e-^/'^^. (16) 

^ ^ TsT \ 2T Us 2TJ 37r2 ^ ^ 

The exponential law, exp(— T/3i?c), reflects the decay of the fiux-dependence of the spin 
fluctuations, while the other factors are due to the temperature dependence of the impurity 
susceptibility. At low temperature (T ^ Ec), and weak magnetic field {ug ^ T), (^2) is 
thus proportional to the inverse temperature, 

(^2) = -,-^E^. (17) 

TT"' Ts I 

The 1/T-divergence is cut off for T <^ Ug, where we find 

B. Strong spin-flip scattering — high concentration of impurities 

It is well known that the cooperon is sensitive to magnetic scattering. As a consequence, 
Eqs. (0) and (0) are only applicable in the limit of weak spin-flip scattering, i.e. for a low 
concentration of magnetic impurities. However, simply inserting the cooperon including 
spin-flip scattering leads to incorrect results, since when calculating the thermodynamic 
potential within the diagrammatic approach, special care has to be taken with symmetry 
factors. 

A typical diagram for the higher order corrections is shown in Fig. ^ The curly line 
represents scattering at magnetic impurities. We consider Ug = 0. In this case, the curly line 
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does not transfer energy, since the thermal average of two spin-operators is time independent, 
e.g. (S'^(r) 5*^(0)) = S{S + l)/3. The summation over all contributions of the type shown in 
Fig. ^ leads to 

,,. 47r;^,n [\uj\ + Dq^ + l/rj 



{n) 



Z — / A^ Z — / ^ 

oo 



1/37 



t^^TT ^^n\ \uj\ + Dq^ + 1/ts + 4/3r,„^ 



(19) 



Performing the n-summation, we find 

leu I + Dq'^ + 2/ts 



(fi) 



2t^27r 



\uj\ +Dq'^ + 2/3ts + 4:/3ts. 



(20) 



|cj| + Dq^ + 1/tJ \ \iu\ + Dq^ + 1/r, + 4/3r,o / 

An expression of this form has been found for the persistent current fluctuations 

and we use results given there in order to evaluate the Fourier expansion of the average 

persistent current. We find 



T $: M in(|,.| + Dq' + 7) = -^^c E ^ cos(47rm0/0o) 

q,uj>0 ^^ """ m=l "^ 



(21) 



for low temperature (T <^ i^c, 7), 



C(7) = e-V^"^'/^= 1 + Vt^V^c + 



7m 

3^: 



(22) 



Using these relations, we can expand Eq. 



as follows: 



in) = ~^^E!f:^cosiAnmH^o), 



T 



n^ 



m=l 



m^ 



(23) 



with 



hr 



C(2/r,) + 3C(2/3r, + 4/3r,o) - C(l/^s) - 3^(1/^. + ^/^rso). 



(24) 



Flux-independent terms have been dropped again. Note that in Eq. (^) the prefactor 
1/T appears which is due to the impurity susceptibility, resulting in rather large Fourier 
components at low temperature. Expanding Eq. (^) in the limit of weak spin-flip scattering, 
but strong spin-orbit scattering, we find hi = —1/6tsEc, recovering the result given in Eq. 
(p!7|) for the grand potential. For vanishing spin-orbit scattering and weak spin-flip scattering. 



the leading term is quadratic in the spin flip scattering rate: hi = {1/tsEc)'^/18. This agrees 
with our result that the contribution to {Q2) which is linear in l/r^ vanishes in the limit 
Us -^ 0. 

For further illustration, we plot the second harmonic of the grand potential as a function 
of the spin-flip scattering rate in Fig. ^ The absolute values of (^2) are comparable in size 
for weak and strong spin-orbit scattering. Over a large range of the spin-flip scattering rate, 
we roughly have 1(^2)1 ~ O.IE'^/T. Without spin-orbit scattering, the persistent current is 
diamagnetic for small magnetic flux, whereas the persistent current is paramagnetic in the 
case of strong spin-orbit scattering. 

In the limit of very strong spin-flip scattering, the persistent current approaches zero. 
Without spin-orbit scattering, the dominant term in Eq. p^) is h^ ~ 3h^{2/3Ts). In the 



limit of strong spin-orbit scattering, hm ~ ~^m(V''"s) is dominant. Thus the persistent 
current, as a function of the spin- flip scattering rate, approaches zero much faster for strong 
spin-orbit scattering than without spin-orbit scattering, as is apparent in the figure. 

C. Strong spin-flip scattering — the Kondo eff^ect 

Perturbation theory in the coupling constant J, as in Eqs. @ and ([T0|), only works 
well far above the Kondo temperature, which is given by Tk ~ ei;'exp(— 1/2| JIA/q). In 
this section, we describe the magnetic impurities in the framework of the single impurity 
Anderson model, and calculate the grand potential to first order in the on-site interaction U. 
In averaging over impurities, we again only keep diagrams with one cooperon. Following the 



ideas of Hewson's renormalized perturbation expansion ||3^, we introduce the renormalized 
Green's functions and a renormalized interaction. This allows calculation of the persistent 
current for temperatures far below the Kondo temperature. We consider low concentration 
of paramagnetic impurities, and Ug = l/r^o = 0. 

The contribution to the grand potential to first order in the Coulomb repulsion between 
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d-electrons is given by 

Und^n^^ = UT^ J2 Gd^{iei)Gdi{ie2). (25) 

Here, Gd is the Green's function of rf-electrons of the non-interacting Anderson model. The 
(i-electrons are coupled to conduction electrons on a disordered ring. The d-Green's function 
averaged over disorder in the conduction band is given by 



Gda{i^) 



d^k '-' 



/ci h 
-——^G^{iuj,k) 



(26) 
(27) 



iu - ed + jJ^ + i6sgn{uj) ' 
where V is the hybridization of the d-electrons with the conduction band, and S the broad- 
ening of the (i-level due to hybridization. Note that the averaged d-Green's function in this 
approximation does not depend on disorder. However the average of the product of two 
Green's functions also involves fluctuations; in averaging Eq. (^Sj) , we keep the diagrams 
shown in Fig. ||, with the result 

M0)) = ^^^E(G'.T(ie))'(G'.i(icc'-ie)))^E(2^Aror)^C+_+_(g,^,e). (28) 

e,aj q 

The fiux-dependence of this equation is due to the low energy singularity of the cooperon. 
Flux independent contributions to the grand potential have been dropped. In order to calcu- 
late the flux-dependence of Q for temperatures below the Kondo temperature, we generalize 
this equation replacing G'^o-(ie) with the exact Green's function Gdai^^n) of the clean Ander- 
son model, and replacing U with the exact irreducible, antisymmetrized four-point vertex 
function. This leads to 

i^m = ^E(^.T(ie))'(G'di(i^-ie))' 

xr}{(ie, ie, icu - ie, icu - ie) E(2vr7Vor)2(:7+_+_(g, u, e). (29) 

In the limit where the Kondo temperature is large compared to the Thouless energy and 
temperature, we can evaluate Eq. ( PP| ) explicitly, since for low energy and temperature, 
analytic expressions for the (i-Green's function and the vertex are known |]32| , |33| : 
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Ti 
with 



G^iie) ^ : — -^ (30) 

le - e^ + yU + 10 

rj}(0,0,0,0)^?7 (31) 



ia = ^i, S = ATK/nw = nAfoV^z, U = ti~5/z^] (32) 

w ~ 0.41 is the Wilson number. The flux-dependence of Eq. (PD|) arises from energies of the 



order of the Thouless energy, so for Tk ^ -Ec, T, we approximate the tZ-Green's functions 
and the vertex by the constant values which they assume for e -^ 0. The flux-dependent 
contribution to the grand potential is then given by 

This is of the same structure as the flux- dependence in the weak coupling limit, see Eqs. 
(10) and (15); thus (f^(0)) contains the following factors: 

1. The scattering amplitude of electrons with the impurities. For T <^ T^ the scattering 
amplitude is given by the unitarity limit for s-wave scattering. The explicit result for 
the scattering rate is l/r^ = 2NsItiNqV. 

2. The magnetic susceptibility of the impurity, which is proportional to the inverse Kondo 
temperature, x = I^'bW /Tk- 

3. A factor which describes the diffusive motion of the electrons around the ring. 
From Eq. (0), we determine the Fourier components of the grand potential (T -^ 0) 

{^2^) = N.-^—,, (34) 

Ti iK 'm 

where we assumed a finite number of impurities; recall that A = l/2A/oV. 

In the presence of spin-orbit scattering we determine the flux-dependence of the grand 

potential from a product of spin-dependent vertex function and cooperon of the form 

o Z^ r^J^^Cai<j2fT3'^4 = y (C+-+- + C-+-+ — C_++_ — C+ — +) . (35) 
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In this combination of the components of the cooperon, which is the singlet component, 
spin-orbit scattering drops out. 

IV. DISCUSSION 

We showed that magnetic impurities contribute significantly to the average persistent 
current in mesoscopic rings. The persistent current is proportional to the impurity suscepti- 
bility. Assuming a Curie law, the persistent current thus becomes large for low temperature 
and dominates the Coulomb interaction contribution to the persistent current. The persis- 
tent current induced by magnetic impurities is large for a wide range of parameters, i.e. in 
the presence or the absence of spin-orbit scattering and also for rather high concentrations 
of impurities. 

Even if the impurity spins are screened due to the Kondo effect, a contribution to the 
persistent current remains. For example for 10'* impurities with T^ ~ IK we find from Eq. 
(p^) (^2) ~ 0.5Ec, when we insert A = 0.2mK which has been reported in the experiment [|I|. 
This is close to the experimental result, where (^^2) ~ O.SEc- However, we have to mention 
that for large concentration of magnetic impurities there are deviations from the linear 
dependence of the persistent current on the concentration of impurities which is assumed in 
Eq. (H. 

A precise calculation of the persistent current as a function of concentration of magnetic 
impurities is not easy, since one has to take impurity- impurity interactions into account. 
Especially, one needs the full energy dependence of the d-Green's functions and of the vertex 
r. In order to estimate the current, we approximate Gd and F by constants, up to a cut-off 
energy of the order of T^-. Thus the dimensionless effective electron-electron interaction due 
to magnetic impurities is given by 

^ = KNj-^\Gd{iuj = 0)\^T][{(},0,0,0) (36) 

^^ Aw w ._. 
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with Hg = Ng/V the concentration of impurities, and Eq. (|^) reads 



UJ 



Including impurity self-energy diagrams in the calculation of the cooperon leads to a renor- 



malization of the diffusion pole |31 



u + Dq'^ ^{1 + v)uj + Dq^, (39) 

and finally, higher order terms renormalize v logarithmically ||28|j3^JT5[| , 



''^''* l + [z//(l + z/)]lnT^/E,- ^^°^ 

However, in any real metal there is the Coulomb interaction and the electron-phonon inter- 
action which complicate our situation further. For low concentration of magnetic impurities, 
the contribution to the persistent current due to Coulomb interaction and magnetic impu- 
rities simply adds, but at a high concentration, this certainly is not the case. In order to 
estimate the current, we approximate the interactions in the conduction band introducing 
a second interaction constant /i. If the Coulomb interaction is dominant, yU is positive and 
the cut-off for the interaction is of the order of the Fermi energy, ujc ^ ep- If the phonons 
dominate, /i is negative with a cut-off of the order of the Debye energy, ojc ~ ojr,. We 
assume that the relation ujc ^ Tk ^ Ec holds. Having two different interactions, we first 
renormalize ^ down to the Kondo temperature, 

f^-^ fi* = — 7—, JT^TT^ (41) 

and then the sum of both is scaled down to E^, with the final result 

in) = 2 ""-^ T V - (421 

^ ^ l + [(^ + ^*)/(l + z.)]ln(T^/E,) ^^^,{l + u)uj + Dq^- ^ ^ 

The Fourier components of the grand potential (and of the persistent current) as a function 
of concentration of magnetic impurities are then given by (T = 0) 

(o ) = "-h^ 4E^ 1 

^ ""^ l + [(z/ + /i*)/(l + z/)]ln(T,^/Ee)7rm3(l + z/)2- ^ ^ 
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With the parameters used above, v ^ 0.5, and there are already strong deviations from a 



hnear dependence of {^2m) on v as given in Eq. (0). 

In conclusion we have shown that magnetic impurities contribute significantly to the 
persistent current. Thus we predict that future experimental studies of persistent currents 
in the presence of magnetic impurities promise to yield most interesting results - and in 
particular a clear test of the theoretical concepts. 
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APPENDIX A: COOPERON IN THE PRESENCE OF SPIN-EFFECTS 

The spin dependence of the cooperon has been discussed at various places in the hterature 
||35| , |28|j3^ . We assume that the rate for scattering at non-magnetic impurities, denoted here 
for clarity by I/tq, is large compared to the spin-flip scattering rate, the spin-orbit scattering 
rate, and the Zeeman splitting, i.e. I/tq ^ 1/ts, \/tso,uJs. The components of the cooperon 
are determined from the equation 



see Fig. [7|; summation over ^ and v is implied here. The bare cooperon is given by 

^%^5 = —^a^^pS + o ( — - — ) ^a7^/3<5- (A2) 

MD 'J \is 'so J 

The quantity 11^,^ is determined from the product over two Green's functions, integrated 
over momentum, 

where /x, z^ = ±1 are spin indices. In the long-wavelength, low-frequency limit the result is 



/ci n 
— -^G^(ie,k)G,(ie-i^,-k + q), (A3) 



n^^ = 27rA/'or [l - r (|cj| + Dcf- + i(/i - z/)wsSgn(cj)/2jJ , (A4) 

where 1/r = I/tq + l/r^ + I/tsq. For the summery of the results we define 

A/'o = |cu| + Dq^ + 2/r, (A5) 

Ari = |cj|+Dg2 + 2/3r, + 4/3r,„, (A6) 

to express the non-zero components of the cooperon in a concise way: 

27rAror2C++++ = 1/iVi (A7) 

27r7Vor2(:7+_+_ = \N^ + A^i - 2ia;,sgn(a;)]/2P (A8) 

27rAror2C+__+ = (A^o - ^xM^V (A9) 
with V = NqNi + ul, and 

c7_+_+ = C7_|___|__, c7_++_ = c7+ — ^, C = C7^_|_^_|_. (AlO) 
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Figure Captions 

Figure [1|: One-cooperon contribution to the fluctuations in the local spin density. 

Figure |^: Second harmonic of the grand potential normalized by the Thouless energy as a 
function of temperature and Zeeman splitting. The spin-flip scattering rate is chosen 
equal to the Thouless energy; no spin-orbit scattering. 

Figure |^: Second harmonic of the grand potential as a function of temperature and Zeeman 
splitting. The spin-flip scattering rate and the Zeeman energy are the same as in Fig. 
0, but here we consider the limit of strong spin-orbit scattering. For comparison we 



also plot (^^2) due to the Coulomb interaction [|T^,|1^. Apparently, the impurity-spin 



induced contribution is dominant at low temperature. 

Figure |^: Higher order Feynman diagram. The curly line represents scattering at magnetic 
impurities. 

Figure |^: Second harmonic of the grand potential as a function of the spin-flip scattering 
rate. 

Figure |^: Graphical representation of the phase sensitive contribution to the grand poten- 
tial in the Anderson model. The double line represents the c?-electron Green's function, 
the wavy line is the on-site interaction U, and the cross denotes the hybridization V. 

Figure |^: Graphical representation of the equation deflning the cooperon. Here, the dashed 
line represents scattering at both magnetic and non-magnetic impurities. 
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